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Abstract. We analyze numerical stability of a recursive computation scheme 
of present value (PV) amd show that the absolute error increases exponentially 
for positive discount rates. We show that reversing the direction of calculations 
in the recurrence equation yields a robust PV computation routine. 



1. Introduction 

The concept of time value of money is so important in finance (see chapter 7 
in [2]), that it is not surprising that the present value (PV) computations are very 
common in finance and accounting. This work focuses on one of the definitions of 
PV, which is similar to equation (7-7) in chapter 7 of [2]. 



Definition 1. Present Value (PV) of future cashflows is defined as follows: 

N r 
m— 1 ^ ' 

where C m is the cashflow amount in month m, N is the month number of last 
non-zero cashflow, r is a monthly discount rate, or required rate of return. Cash 
outflows are negative, and cash inflows are positive. 

Definition 2. Net Present Value (NPV) of cashflows is defined as follows : 
(1.2) NPV{r) = C + PV(r) 

where an initial cash outflow is Co, i.e. month zero. Cash outflows are negative, 
and cash inflows are positive. 

Of course, the most trivial method of calculation of PV is by direct summation 
of discounted cashflows n!ffim i i- e - J us ^ ^ e m ^ ne definition. However, it is easy 
to to derive a recurrence equation. First, let us define a present value of the subset 
of a given cashflow as follows. 



Definition 3. The partial PV of a cashflow is defined by equation: 

N 

(1-3) PV k = £ 

, , , (1 + r) 



(1 + r) m 

m=fe+l v ; 

Obviously, PV$ is PV of an entire cashflow excluding initial payment, i.e. PVq = 
PV(r). From the time value point of view, PVk is nothing but the present value of 
the remaining cashflows at month k (exclusive) . Knowing a partial P V in previous 
month, it is easy to calculate it for the current month using the recurrence equation: 

(1.4) PV k = (1 + r) ■ PV k ^ - C k 



Date: November 10, 2006. 

1991 Mathematics Subject Classification. Primary 65G50, Secondary 65-04 . 
Key words and phrases, present value pv recursive recurrent numerical instability. 

1 



2 



ARGYN KUKETAYEV 



This formula can easily be derived as follows: 

n n 

m=fc+l V ' m=k y ' 

2. Recurrent PV computation use case 

We believe that the equation 11.41 is used in practice. Consider the Table 1 in 
Appendix B of lj . The document was published by Financial Accounting Standards 



Board (http:/ /www. fasb.org/) and it sets a guideline on amortization of certain fees 



and costs associated with loan origination. We refer to chapter 15 and page 688 in 
[2] for discussion of amortization, amortization schedule and loan origination fees. 
Table [T] is essentially a fragment of the above mentioned table with insignificant 
modifications. We used Google Spreadsheet application [3] to create this table by 
implementing FAS 91 rules in the worksheet. We show it here in order to make a 
point that the recursive PV computation algorithm is actually used in practice and 
also to explain why it is needed. 

Table 1 . Amortization Based on Contractual Payment Term 





(1) 


(2) 


(3) 


(4) 


(7) 


Year 


Cash (Out) 


Stated 


Amortization 


Interest 


Carrying 




Inflow 


Interest 




Income 


Amount 




-$98,000.00 








98,000.00 


1 


$16,274.54 


10,000.00 


263.52 


10,263.52 


91,988.98 


2 


$16,274.54 


9,372.55 


261.44 


9,633.99 


85,348.43 


3 


$16,274.54 


8,682.35 


256.18 


8,938.53 


78,012.42 


4 


$16,274.54 


7,923.13 


247.10 


8,170.23 


69,908.11 


5 


$16,274.54 


7,087.99 


233.48 


7,321.46 


60,955.04 


6 


$16,274.54 


6,169.33 


214.48 


6,383.81 


51,064.31 


7 


$16,274.54 


5,158.81 


189.15 


5,347.96 


40,137.72 


8 


$16,274.54 


4,047.24 


156.38 


4,203.62 


28,066.80 


9 


$16,274.54 


2,824.51 


114.92 


2,939.43 


14,731.69 


10 


$16,274.54 


1,479.50 


63.34 


1,542.85 


0.00 




Total 




$2,000.00 








amortization 











We have to introduce the concept of the internal rate of return (IRR) (see page 
202 in |2J) in order to deconstruct Table [TJ 

Definition 4. Internal rate of return (IRR) of cashflows is a rate r, which makes 
NPV(r) equal to zero and is defined by the following equation: 

N 



(2.1) NPViirr) = V - — "I" . = 

y ' ^— ^ (l + irr) m 

In the Table Q] column (1) contains cashflow amounts for each year. In FAS 
91 examples years are used instead of months. It does not change anything in 
principle if we assume that years are months, for the time being. In our notation 
C* = -98, 000.00 and N = 10. Next, the IRR is calculated for this cashflow and it 
is approximately 10.47%. Now, column (7) (carrying amount) is actually PVk(irr), 
e.g. carrying amount in year 9: 

NPV^rr) = -^2_ = 16 ' 27454 = 14,731.69 
y 1 l + irr 1 + 0.1047 

The following is an excerpt from the description to Table 1 in PQ: 



C r , 
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Computations: 

Column (1): Contractual payments 

Column (2): Column (5) for prior year x the loan's stated interest 
rate (10%) 

Column (3): Column (4) - Column (2) 

Column (4) : Column (7) for prior year x the effective interest rate 
(10.4736%) & 

Column (5): Column (5) for prior year - (Column (1) - Column 
(2)) 

Column (6): Initial net fees - amortization to date 
Column (7): Column (5) - Column (6) 

In this description IRR is referred to as effective interest rate. Also, in our copy 
of this table we omitted columns (5) and (6) in order to save the space. We leave 
it to the reader to prove that carrying amount is equivalent to PV (definition 1 1.3|) . 
and that it is computed recursively (equation ll.4|) starting from the top of the table. 
We shall call this algorithm forward recursive PV computation scheme, as opposed 
to backwards scheme, which will be described later. The table itself is often called 
amortization schedule. 

Why not compute carrying amounts directly using PV formula 11.11 / Consider 
the following use case. A company may have a huge portfolio of loans, maybe 
500,000 loans. Long-term loans, such as mortgages, may have scheduled terms 
of 30 years (360 months) or even more. Often more than one fee is amortized. 
Direct calculation of NPV for each month in the schedule yields excessive redundant 
summation of the same discounted cashflow amounts over and over. Therefore, it 
is very desirable to have a recursive PV computation method at hand. 

3. Numerical instability of forward recursive PV computation for 

POSITIVE IRR 

Let us analyze absolute error accumulation when using forward PV computa- 
tion(eq. 11.40 . The following equations describe accumulation of rounding errors in 
floating-point addition/subtraction and multiplication. 

(3.1) A affi6 ~A a + A 6 

(3.2) A a ® b ^ A a ■ \b\ + A b ■ \a\ 

These equations are commonly used in analysis of error propagation. Essen- 
tially, they are equivalent to the notion of accumulation of absolute errors in sum- 
mation/subtraction operations, and accumulation of relative errors in multiplica- 
tion/division operations. One can derive them using partial derivatives as follows 
(see 2.11 in g]) 

d d 

(3.3) A a ® b « |— (a • 6)| • A« + |^(a • 6)| • A 6 

Applying these equations to equation 1 1.41 we can deduct the following 

(3.4) A PVk ^\l + r\-A PVk _ 1+ \PV k ^\-A r + A Ck 

We shall assume that Ac k = 0, i.e. cashflow amounts are precisely accurate. 
Also, PVo = Cq, because we used irr as a discount rate. Therefore, Apy = and 
the following is true 

(3.5) A PVl «PV - Ar = | Co | • A r 

(3.6) A PV2 » |1 + r\ ■ \C \ ■ A r + \PV k -i\ ■ A r 
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(3.7) A PV3 w |1 + r\ 2 • | Co | • A r + |1 + r\ ■ \PV k -i\ ■ A r + \PV k -i\ ■ A r 

and so on. It is easy to show by simple mathematical induction that the absolute 
error of present value has a term, which increases exponentially with k for positive 
rates r: 

(3.8) A PVk « jl + rf- 1 • |C | ' A r 

Note that A r can not be zero, because it is computed by numerically solving equa- 
tion [23] Therefore, we state that the forward recursive PV computation is numer- 
ically unstable for r > 0. 

In order to demonstrate our finding, let us compute PVk for every year in Table[T] 
using direct summation formula ll.3l then compute the rate of year-on-year increase 
of the difference between this value and the carrying amount. The results are shown 
in Table [J]. If we assume that direct PVk computation by equation 11.31 is a better 
estimate of true present value, then column (9) confirms the correctness of our error 
analysis as per equation 13.81 



Table 2. Amortization Error Accumulation Exhibit 





(8) 


(9) 


Year 


PV k (as in OK 


% increase from 




- Column (7) 


previous year 


1 


0.000000000509317 




2 


0.000000000567525 


11.43% 


3 


0.000000000640284 


12.82% 


4 


0.000000000698492 


9.09% 


5 


0.000000000778527 


11.46% 


6 


0.000000000858563 


10.28% 


7 


0.000000000945874 


10.17% 


8 


0.0000000010441 


10.38% 


9 


0.000000001149601 


10.10% 


10 


0.000000001269882 


10.46% 



4. Backward recursive PV computation is robust for positive IRR 

Let us look at equation [Lj] and its absolute error formula [3. 81 There seem to be 
two problems. 

• Every next present value is being multiplied by factor 1 + r, so the absolute 
error of PVk should increase every month by percentage equal to r for 
positive rates. 

• PVk itself is decreasing toward the end of the schedule until it finally reaches 
zero. 

Combined, these two factors should magnify the relative error of the present at 
every step of the calculations. 

A heuristic suggests considering to reverse the recurrence equation 11.41 

(4.1) PVk-i = ^— (PV k + C k ) 

1 + r 

We shall call the new recurrence equation a backward recursive computation 
method. It seems to address both mentioned issues. 



^Note that the discount rate in this sample is 10.47%. 
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• Its next step is evaluated by division of the previous value by 1 + r, subse- 
quently the absolute error of the present values should not be magnified. 

• The sum in the second factor generally increases on every step, i.e. unlike 
11.41 it does not add to relative errors of computed values. 

Like we did it earlier for forward method, let us conduct a simple error prop- 
agation analysis. First, we shall define the absolute error of using equation 



(4-2) A i w r — 

17 — (l + r )2 

Now, we can estimate the error of the backward recursive computation as follows, 
again, assuming Ac fc = 

( 4 - 3 ) A pv^ ~ + (i + r)2 

Note, that PVn = and Apy N = 0, because we are using IRR as a discount 
rate. Now we can find out absolute errors of first steps of backward computation. 

(4-4) Ap^ « ° + = 



(4.5) A 



pv N - 



\l+r\ (l + r) 2 
|PW_ 1 |-A r _ |PVw_i|-A r 



|l + r| (l + r) 2 (l + r) 2 



(Art A \PVN-l\-Ar , \PV N -2\-A r 

(4-6) A Pyjv _3 ~ | 1 + r |3 + (i + r)2 

and so on. We can generalize these equations for positive discount rates r using 
mathematical induction as follows 

(4.7) A PVk J-^±^ + 0^ PV ^-^ 



(l + r) 2 V |l + r| 3 

As expected, for positive discount rates r, there are no exponentially increasing 
terms. The only increasing factor of Apy k is of order PVk+i- However, it does not 
magnify the relative error, because we know that generally PVk > PVk+i- 

For negative discount rates r, Apy N _ k will have terms like ][PV ^^ll n ^ r ■ Similar 
to forward computation with positive r, the factor | 1+r |fc-m increases exponentially 
with k. However, as it was stated earlier, PVk decreases toward the end of an 
amortization schedule. Therefore, in backward computation PV/v-m+i in numera- 
tor should compensate, at least partially, the shrinking denominator |1 + r \ k ~ m . 

In order to test our findings, we produced tables [T] and [3 using backward com- 
putation scheme. As expected, we did not observe exponential increase of absolute 
errors in recursively computed present values. In order to deal with this instability, 
as it is suggested by theoretical error propagation analysis, it is necessary to in- 
crease an accuracy of IRR calculation, i.e., in other words, decrease A r . However, 
the exponential character of absolute error accumulation and the finite precision of 
floating-point arithmetics limits the remedial power of this approach. 

Conclusion 1. We demonstrated numerical instability of forward recursive present 
value computation method. We showed that this method is well known and might be 
used in accounting and financial software. We stated that backward recursive present 
value computation method is robust, and suggested to consider using it instead of 
forward method for positive discount rates. 
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